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Abstract. We describe algorithms, based on Avis and Fukuda's reverse search paradigm, for listing 
all maximal independent sets in a sparse graph in polynomial time and delay per output. For bounded 
degree graphs, our algorithms take constant time per set generated; for minor-closed graph families, 
the time is 0{n) per set, and for more general sparse graph families we achieve subquadratic time 
per set. We also describe new data structures for maintaining a dynamic vertex set S in a sparse or 
minor-closed graph family, and querying the number of vertices not dominated by S; for minor-closed 
graph families the time per update is constant, while it is sublinear for any sparse graph family. We 
can also maintain a dynamic vertex set in an arbitrary m-edge graph and test the independence of 
the maintained set in time O(v'm) per update. We use the domination data structures as part of our 
enumeration algorithms. 

1 Introduction 

In this paper we revisit the classical combinatorial enumeration problem of generating all maximal 
independent sets of a graph, or equivalently generating all cliques in the complement of the graph. 
Many algorithms for this problem are known; see, e.g., [9-13,15,17,18,21,22], or the survey by 
Bomze et al. [3]. However, all take at least quadratic time per generated independent set even 
on sparse graphs. There has also been work on maximal independent set enumeration with time 
bounds proportional to the maximum possible number of such sets [4,8,20], but these algorithms 
may take a large time per set on instances with few maximal independent sets. We find more 
efficient maximal independent set generation algorithms for many graph classes including bounded 
degree graphs, minor-closed graph families, and subgraph-closed sparse graph families. Although 
graphs in these classes have polynomially many cliques [5] , generating all maximal independent sets 
is more difficult for these graphs as there may be exponentially many such sets. 

Our maximal independent set algorithms require us to quickly test whether certain subsets of the 
vertices of the graph dominate the remaining vertices, so we also consider the problem of dynamic 
dominance testing. We provide very efficient algorithms for this problem on minor-closed graph 
families, taking constant time per test. We also find sublinear algorithms for dynamic dominance 
in more general sparse graph families. As a simple consequence of the latter result we show how to 
test independence of dynamic vertex sets in arbitrary graphs. 



2 Reverse Search 

Reverse search is a powerful paradigm for enumeration developed by Avis and Fukuda [1,2] and 
applied by them to many enumeration problems including listing all vertices of a convex polytope, 
triangulations of a planar point set, vertices or cells of a hyperplane arrangement, spanning trees in 
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graphs and non-crossing spanning trees in the plane, connected induced subgraphs of graphs, and 
topological orderings of directed acyclic graphs. 

To use reverse search, one needs a parent operation on the objects being enumerated, such that 
any two objects can be transformed to a common canonical object by repeated parent operations. For 
instance, for triangulations, the parent operation finds two triangles sharing an edge and maximizing 
the sum of the two angles opposite the shared edge, and replaces them by two other triangles with 
the same union; repeated application of this flip operation reduces any triangulation to a canonical 
triangulation, the Delaunay triangulation of the input. For convex polytopes, the parent operation 
can be any simplex method pivot rule, for an appropriate linear objective function, and the canonical 
object is the vertex minimizing the objective function's value. We form a rooted tree, having as its 
vertices the objects to be enumerated, with edges between each object and its parent, and having 
the canonical object as root. Reverse search is simply depth first tree traversal applied to this 
tree. In order to perform this traversal, one must be able to identify the children of each node in 
the tree. Typically, a polynomially-sized superset of the children is identified (e.g. in the case of 
triangulations, this superset could be all triangulations reachable from the given one by a single 
flip). We can then test whether any member of the superset is a child of the current object by 
applying the parent operation to it and comparing the result to the current object. 

Thus, reverse search performs the following steps for each object: identify a superset of the 
object's children, apply the parent operation to test each member of the superset, and search 
recursively each child found by this test. In this basic form, the time per object is the product of 
the number of items in the superset with the time to generate each member of the superset and 
apply the parent operation to it. As we shall see for our maximal independent set enumeration 
algorithm, however, reverse search may be sped up in various ways, for instance by maintaining 
data structures that allow us to perform the parent tests more quickly. 

3 Basic Independent Set Generation Algorithm 

To apply reverse search to maximal independent set generation, we need an appropriate parent for 
each maximal independent set S. Assume we have ordered the vertices into a sequence; this allows 
us to compute the lexicographically first maximal independent superset of any independent set by 
considering the vertices one at a time, in order, and adding each vertex to the set whenever it is 
independent of the vertices already in the set. The lexicographically first maximal independent set 
(LFMIS for short) is the lexicographically first maximal independent superset of the empty set. 
To find the parent of S, let v be the vertex of LFMIS \ S occurring earliest in the sequence, let 
N = N{v) n S, and let parent(5) be the lexicographically first maximal independent superset of 
{S U {v}) \ N. Each parent operation increases the length of the initial prefix of shared vertices 
between S and LFMIS, so after at most |LFMIS| repetitions of the parent operation, the canonical 
maximal independent set LFMIS will be reached. Therefore, this parent operation satisfies the 
requirements of the reverse search paradigm. 

For any vertex v in LFMIS, let later(f ) be the set of vertices in V{G) \ LFMIS for which v 
is the adjacent LFMIS vertex that appears earliest in the ordering. The sets later (w) partition 
V{G) \ LFMIS. If S is maximal independent, and parent(S') is the lexicographically first maximal 
independent superset of {S U {v}) \ N, then N C later(f), for no vertex in N could be adjacent 
to an earlier vertex in LFMIS than v without violating the assumed independence of S. Thus, to 
search for the children of a maximal independent set S, we need only consider sets N C later(f ) 
for some v in the initial prefix of shared vertices between S and LFMIS. 

The recursive search pseudocode in Listing^ with an initial call to search (LFMIS), instantiates 
the reverse search paradigm for the maximal independent set problem. The collection of potential 
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def search(S): 
output S 

for each vertex v in the ordered sequence: 
if V is not in LFMIS: 

continue 
if V is not in S: 

break 

for each nonempty independent subset N of later(v): 
T = (S union N) \ (neighbors of N) 
if T is maximal and parent(T) == S: 
search (T) 

Listingl. Recursive reverse search for maximal independent sets. 



children generated by each call to search may include sets that are not maximal independent, so 
we must check that each generated set T is maximal as well as that its parent is S] this problem 
of testing maximality is a large part of the difficulty in making our reverse search algorithms 
efficient. To reduce the space occupied by the call stack, we generally prefer a version of the search 
procedure in which the recursive calls have been unfolded into the procedure itself. This unfolding 
can be done without need for an auxiliary stack, at the expense of some code complexity and an 
additional parent computation per generated set, as shown in Listing |21 

4 Minor- Closed Domination 

To make our search procedure efficient, we need a fast way of handling its maximality tests. The 
sets that must be tested are automatically independent, and an independent set is maximal if and 
only if it dominates the graph; that is, all vertices are either in or adjacent to the set. Therefore, 
we make use of a more general data structure for testing domination. We describe here such a data 
structure for minor-closed graph families. That is, given graph G from a minor-closed graph family 
J^, we wish to maintain a set 5 C V{G)^ subset to insertions and deletions of elements of 5, and 
answer queries requesting the number of vertices in G that are not dominated by S. We show how 
to do this in constant time per update or query, and linear space and preprocessing time. We will 
use this data structure later as part of our maximal independent set enumeration algorithm. 

Lemma 1. For any minor closed graph family T , graph G ^ T , and nonempty Q C V{G), define 
an equivalence relation ~ on V{G) \Q by u ^ v if and only if N{u) H Q = N{v) n Q. Then the 
number of equivalence classes of ^ is 0{\Q\), where the constant of proportionality depends on T 
but not on G. 

Proof. Suppose for a contradiction that we can find G £ J- and Q C V{G) where the number of 
equivalence classes of ~ is an arbitrarily large multiple of \Q\. For any such Q and G, form a minor 
Gq of G by the following process: let the vertices of Gq be the same as the vertices of Q. Consider 
(sequentially, in arbitrary order) the equivalence classes of ~, choose a representative vertex v for 
each equivalence class, and form an edge in Gq by contracting a path through v between two vertices 
in N{v) n Q, unless all such paths connect pairs of vertices that are already connected in Gq. At 
each step of the construction process outlined above, Gq is a minor of G; therefore, it belongs to J^, 
and using the known sparsity of minor-closed graph families [14] we can show that there are 0(|Q|) 
cliques in the partially constructed minor Gq [5]. An equivalence class represented by v can only 
fail to add an edge to Gq iif N{v) n Q already forms a clique, so only 0(|Q|) equivalence classes 
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def search(S): 

Iseq = subsequence of vertices v with nonempty later (v) 
while True: 
output S 

V — first vertex in Iseq 

N = first nonempty independent subset of later(v) 
while True: 

T = (S union N) \ (neigiibors of N) 
if T is maximal and parent(T) == S: 
S = T # unfolded call to search (T) 
break 
while True: 

if N != last nonempty independent subset of later(v): 

N = next nonempty independent subset of later (v) 

break 
if v != last in Iseq : 

V — next in Iseq 

if V in S: 

N — first nonempty independent subset of later (v) 
break 
if S == LFMIS: 

return 
P = parent(S) 

V = first vertex from Iseq in P \ S 
N = S \ P 

S = P # unfolded return to search (P) 

Listing2. Nonrecursive reverse search for maximal independent sets. 



fail to contribute an edge. Thus, if the number of equivalence classes could be an arbitrarily large 
factor times \Q\, we could form arbitrarily dense minors Gq, contradicting the known bounds on 
density of minor-closed graph families [14]. □ 

One can prove more directly that planar graphs have at most max(6|(5| — 9, 2|Q|) classes. 

Our data structure consists of the vertices of G together with a linear number of additional 
supervertices. To form the supervertices, we form a sequence of graphs Gq = G, Gi, G2, . . . , as 
follows. Let A he a constant, depending on graph family but not on the particular graph G. To 
form Gi from Gj-i, let Qi be the set of vertices in Gj_i with degree at least A. Partition Gj_i \ Qi 

into the equivalence classes of ~ and form a single supervertex for each such class. Form the edges 
of Gi from the induced subgraph for Qi, together with an edge from a supervertex s to h £ Qi 
whenever a vertex v in the equivalence class corresponding to s is connected to h in Gj_i. 
We can make the following observations and further definitions about this process. 

— Each Gi is isomorphic to a subgraph of G (formed by choosing a representative vertex for each 
supervertex and omitting edges between pairs of supervertices), so it belongs to J^. 

— By Lemma n |y(Gj)| = 0{\Qi\), and by choosing A sufficiently large we can make \Qi < 
e|y(Gj_i| for any constant e > 0. Thus, we can ensure that each graph in the sequence is 
smaller by a constant factor than the previous one. 

— Each supervertex has degree at most Z\ — 1, so each vertex in Qi is an original vertex of G. 

— The sequence of graphs Gj terminates only when the remaining graph forms a single supervertex 
of degree zero. 

— Each equivalence class of ~ contains at most one supervertex of Gj_i with the same degree, 
item If a supervertex v of Gi contains a supervertex u of Gi_i, and has the same degree, we 
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consider the two to have the same identity, and represent them by the same object in the data 
structure. 

— Define the level of a vertex or supervertex v to be the largest i such that v G V{Gi). Each 
original vertex of G belongs to at most A — 1 supervertices, which (if ordered by level) have a 
decreasing sequence of degrees. 

Our data structure consists of the following information. 

— The graph G and set S to be maintained, and the set of vertices and supervertices in all graphs 
Gi constructed as described above. 

— For each vertex or supervertex v of level i, other than the degree-zero supervertex, a pointer 
sy{v) to the supervertex corresponding to the equivalence class v belongs to in Gj+i. 

— For each vertex or supervertex of level i, a count nadj('y) of the number of adjacent original 
vertices of G that belong to S and have level at most i. 

— For each vertex of G with level i, a list nbr(u) of adjacent vertices and supervertices in Gi. 

— For each vertex or supervertex, a number undom(v). For a vertex of G, undom(f) = 1 when 
V G y{G) \ S and undom(v) = when v e S. For a supervertex, 

undom(f) = ^ undom(i(;). 

{w I sv{w)=v A nadj(«))=0} 

Theorem 1. For any minor- dosed graph family T and graph G & J^, the data structure above 
requires 0{n) space and can be constructed in 0{n) time. We can query the number of undominated 
vertices in G in 0(1) time, and insert or delete vertices in S in 0(1) time per update. All constants 
of proportionality in these bounds depend only on T and not on G. 

Proof. The only nontrivial step in the construction of each graph Gi is finding the equivalence 
classes of ~^ which can be done by bucket sorting in 0(n) time. The times for constructing the 
whole sequence of graphs Gi add in a geometric series to 0(n). 

To query the number of undominated vertices, return undom(z) where z is the degree-zero 
supervertex. We say that supervertex v is reachable from vertex u ii u == u or u is reachable from 
sv(ii); z is reachable from every vertex. If v is undominated, nadi{w) = for each w reachable 
from V, and v contributes one to each reachable supervertex. If v is in S, it contributes zero to 
each reachable supervertex. If v is not in S, but is dominated by a neighbor u £ S, let vj be the 
supervertex containing v at the same level as u; then nadj {w) > and v does not contribute to any 
supervertices reachable from w. Therefore, undom(z) is the number of undominated vertices in G. 

To insert a vertex v to S, increment nadj(«;) for every w in nbr(?;), and (if this causes nadj(';«) to 
change from to nonzero) update undom(x) for every x reachable from w. Also update undom(f ) 
and undom('u) for every v reachable from v. Each update takes constant time and a constant number 
of updates are performed, so the total time is constant. Deletions are handled similarly. □ 

5 Speirse Domination 

We now consider domination data structures for more general sparse graph families. We say that 

graph G is k-orientable if we can orient the edges of G in such a way that each vertex has out-degree 
at most k. Equivalently (by Hall's theorem), this condition asserts that every subgraph H C G has 
at most k\V{H)\ edges. Any subgraph-closed family of graphs with 0(n) edges per n- vertex graph 
is A;-orientable for some constant k; for instance, planar graphs are 3-orientable, and a 3-orientation 
of a planar graph may be found in linear time [6] . 
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Our data structure for domination in sparse graphs resembles that for minor-closed graph 
families, but differs in detail. Given A:-oriented graph G and vertex set Q C V{G), define equivalence 
relation ^ on V{G) by u u if and only if N^{u) r\Q = N^{v) n S", where maps a vertex to 
its outgoing neighbors. If s is an equivalence class of '~ containing f , let N~^{s) = N'^{v) H Q. We 
choose a (nonconstant) value A, let Q be the set of vertices with degree at least A, and create a 
supervertex for each equivalence class of In our data structure we store the following data: 

— The dynamic set S for which we wish to maintain dominance information. 

— For each vertex of G, the supervertex corresponding to its equivalence class. 

— For each vertex v of G, the number lowdom(i;) of edges uv where u belongs to S, and where 
either uv is oriented from u to v oi u has degree less than A. 

— For each supervertex s, the number nundom(s) of vertices v in its equivalence class which are 
not in S and for which lowdom(u) = 0. 

— For each supervertex s, the number hidom(s) = |A^^(s) n S\. 

— The sum of nundom(s), summed over those supervertices for which hidom(s) = 0. 

Theorem 2. For any k-orientable graph G, the data structure above requires 0{n) space and can 
he constructed in polynomial time. We can query the number of undominated vertices in G in 0(1) 
time, and insert or delete vertices in S in 0(n^~^/'^) time per update. 

Proof. The query answer is given by the overall sum of nundom(s). 

To insert vertex v into S, increment the counts lowdom(t') of all outgoing neighbors of v, and 
(if v has low degree) all incoming neighbors of v. Update nundom(s) for each supervertex s that 
has an equivalence class containing v or one of these neighbors. If v has high degree, increment 
hidom(s) for all supervertices s with v E N~^{s). Whenever we change nundom(s) or hidom(s) we 
update the overall sum. The process of deleting a vertex is similar. 

Inserting or deleting a low degree vertex takes time 0{A), and inserting or deleting a high 
degree vertex takes time 0{{n/A)^~^) since there are that many supervertices associated with a 
fixed high degree vertex. By choosing A = n^-^/f' we achieve the stated bounds. □ 

Dominating sets in A;-oriented graphs can be used to model other problems including inde- 
pendent sets, matching, fc-SAT, and constraint satisfaction. As an example we show how to test 
independence in general graphs. 

Corollary 1. We can maintain a set S of vertices in an arbitrary m-edge graph G, and test the 
independence of S, by a dynamic data structure that takes time 0(1) per test and 0{y/m) per 
insertion or deletion in S, and uses linear space and preprocessing time. 

Proof. Form a graph G' having vertices corresponding to sets of 0, 1, or 2 vertices in G; we include 
as 2-vertex sets in G' only the sets of endpoints of edges in G. Connect two vertices in G' by an 
edge whenever the corresponding sets differ by a single element. Then G' can be 2-oriented by 
orienting all edges from larger sets to smaller sets. The subset S is independent in G if and only 
if {0} U {{v} : V ^ S} dominates all vertices of G' , so independence in G can be tested by our 
domination algorithm for the 2-oriented graph G' . □ 

The same result can be achieved more directly: In any fc-oriented graph, we can maintain the 
number of adjacent pairs in a dynamic set S, and therefore determine the independence of S, by 
a simple data structure that stores for each vertex the number of incoming edges from vertices in 
S, in time 0{k) per update. The corollary follows since any m-edge graph can be 0(-v/m)-oriented. 
However, the proof we have given for Corollary ^ provides some evidence that dominance is strictly 
harder to maintain than independence, since any improvement to Theorem|2lfor 2-orientable graphs 
would lead to a corresponding improvement to Corollary ^ for arbitrary graphs. 
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6 Sparse Independent Sets 

A graph G is k-degenerate [16, 19] if its vertices can be ordered in such a way that, for each vertex, 
the number of neighbors occurring later in the ordering is at most k. Equivalently, each subgraph 
of G has a vertex with degree at most k. This parameter is also known as the inductiveness or 
the Szekeres-Wilf number of G. A /c-degenerate ordering of G, if one exists, can be found by a 
simple greedy algorithm in linear time [16]. It is known [14] that all minor- free graph classes have 
bounded degeneracy. A /c-degenerate graph is clearly A;-orientable (simply orient each edge from 
the earlier to the later vertex in a /c-degenerate ordering) and conversely a fc-orientable graph must 
be at most 2/c-degenerate, so graphs of bounded degeneracy are the same as the graphs of bounded 
orientability considered in the previous section. 

Theorem 3. Let G be a k-degenerate graph in which we can maintain a dynamic set S and test 
whether S dominates V{G) in time T{n) per insertion or deletion to S. Then, we can list all maximal 
independent sets in G, in time 0{nT[n)) per generated set and polynomial delay. The space required 
by the algorithm is 0{n) plus a single instance of the dynamic dominance data structure. 

Proof. We use the /c-degenerate ordering as the vertex ordering for our reverse search algorithm. 
Therefore, each set later(t>) will have at most k vertices, each vertex in G will participate in 0(1) 
sets C later (u) throughout a call to search, and the sum of the numbers of neighbors of sets 
A^ will be proportional to the number of edges in the graph, which is 0{n). To save space, we 
unfold the recursive calls of the search into a nonrecursive version of the algorithm, as described in 
Listing 121 The delay bound follows since we can at most test 0{n) potential children each at 0{n) 
levels of the recursion before outputting another set or exiting the search. 

As the algorithm progresses, we maintain a dynamic dominance data structure for the current 
set S, which we use to test each successive set T for maximality with a number of updates propor- 
tional to the size of the set of neighbors of N; thus, throughout a call to the recursive version of the 
search algorithm, the number of data structure updates is 0{n). In the nonrecursive search, we also 
update the data structure whenever we change the set S by an unfolded recursive call or return; 
the total number of data structure updates caused by these changes is again 0(n) per output set. 

Finally, we must consider the time taken to compute parent (T) for each potential child T 
considered by the search algorithm. In these computations, parent (T) is the lexicographic maximal 
independent superset of the independent set (T U {v}) \ N. We also know (from the maximality 
of r) that the only vertices that can be added in forming the lexicographic maximal independent 
superset are neighbors of A^. To perform this computation efficiently, we maintain a simple data 
structure that stores a count in each vertex of the number of incoming edges from vertices in S. 
Each change to S causes k counts to be updated, in time 0{1). Then, we modify this data structure 
to count incoming edges from T instead of S, and use it to compute the lexicographic maximal 
independent superset of T in time proportional only to the number of neighbors of N. Therefore, 
all parent computations can be done in time 0{n) per generated set. □ 

Corollary 2. Let T he a minor-closed graph family, and let G he an n-vertex graph in T . Then 
we can list all maximal independent sets in G, in time 0{n) per set, space 0{n), and polynomial 
delay, where the constants of proportionality depend on T hut not on G. 

Corollary 3. We can list all maximal independent sets in any n-vertex k-oriented graph in time 
0{n'^^'^/^) per set, space 0{n), and polynomial delay. 
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7 Bounded Degree Independent Sets 



We now briefly describe our algorithm for bounded degree graphs. The key observation in this case 
is that, for a given maximal independent set S, vertex v € S D LFMIS, and C later (t;), the 
associated set T is a child of S if and only if the following three conditions hold: 

1. T is maximal. 

2. S is the lexicographically first maximal independent superset of (T U {v}) \ N. 

3. Vertex v is the earliest vertex of the sequence in LFMIS \ T. 

The first two conditions depend only on the inclusion or exclusion in S" of a constant number 
of vertices within distance 0(1) of v. In particular, T is non- maximal if and only if some vertex 
within distance three of v can be added to T, which can be tested by examining all vertices 
within distance four of v. If T is maximal, the computation of the lexicographically first maximal 
independent superset of (T U {v}) \ N can be done by examining only vertices within distance two 
of V. 

We say that a pair (f , A^) is fertile for a set S if the set T generated from S using v and N 
passes the first two of the three conditions listed above. The third condition can be rephrased as 
stating that v belongs to the initial common prefix of LFMIS and S. As our algorithm progresses, 
we maintain the following information: 

— The position of the last vertex in the initial common prefix of LFMIS and S. 

— The set of all fertile pairs {v, N) where v occurs before the last initial common vertex, stored 
as a dictionary mapping each vertex v to the sets A^ that form fertile pairs for it. 

The pairs in the set maintained by the algorithm give exactly the children of the current maximal 
independent set S. Each child differs from S in a constant number of vertices, and so the set of 
fertile pairs for the child also differs by a constant. However, the position of the last vertex can 
differ dramatically between S and its children. In order to keep the changes to the set of fertile 
pairs gradual, we modify our search algorithm so that it processes the children of S in the reverse 
of the vertex ordering on the vertices v. 

Thus, the maximal independent sets in G can be generated by the algorithm described in 
Listing 121 Each time we change S to form one of its children, we can potentially affect fertile pairs 
for vertices within distance four of the change; for each such vertex w occurring no later than 
last_common in the vertex sequence, and each independent subset A^ of later (it;), we test whether 
the change to S has caused {w, N) to start or stop being a fertile pair, and if so add or remove A^ 
to or from fertile_pairs [w]. There are 0(1) pairs {w, N) tested per change to S, so the total time 
per child is 0(1). 

We have not described what order to use for the vertices, because any ordering will work for 
the correctness of the algorithm and its asymptotic analysis. However, a /c-degenerate ordering for 
the minimum possible k may be preferable to other orderings, because it reduces the number of 
subsets of later (u) that need to be considered for each v and thereby reduces the constant factors 
in our analysis. 

One complication with the analysis of the algorithm above is the question of how we maintain 
or sort the vertices v considered by the outer loop, so that they are considered in reverse order. For 
this analysis, we assume a simple comparison sorting algorithm that sorts these vertices in time 
0(A;log A;), where k is the number of vertices to be sorted. 

Lemma 2. // there are k vertices in the fertile_pairs data structure for maximal independent 
set S, then Q{k) of the children of S have Q{k) vertices in their respective fertile_pairs data 
structures. 
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def search(S, last_common, fertile_pairs ): 
output S 

for V in fertile_pairs , in reverse order by vertex sequence: 
later .subsets = fertile_pairs [v] 
delete v from fertile_pairs 
last_common = predecessor of v in LFMIS 
for N in later_subsets : 

S = (S union N) \ (neighbors of N) 

update fertile_pairs from changes to S 

search(S, last .common, fertile_pairs ) 

S = LFMISS((S union {v}) \ N) 

update fertile_pairs from changes to S 

initial call : 
S = LFMIS 

last_common = last vertex in S 
F = fertile pairs for S 
search(S, last_common, F) 

Listings. Reverse search for maximal independent sets in bounded degree graphs. 



Proof. The changes in the maximal independent set from S to its children lead to 0(1) changes to 
the fertile_pairs [w] data structure, and each time we consider a vertex v in the outer loop for S 
we remove only that vertex from the data structure. So, for each of the first k/2 vertices considered, 
there remain k/2 — 0(1) vertices in the data structure for the children. □ 

Theorem 4. Let G be a graph with maximum vertex degree 0(1). Then, we can list all maximal 
independent sets in G, in time 0(1) per generated set, space 0{n), and polynomial delay. 

Proof. We use the recursive version of the algorithm described in Listing |31 which modifies the set 
S and the data structure fertile_pairs [w] in-place and shares the modified structures with each 
recursive call. The space for these structures is 0(n), the other space used by the algorithm per 
call is 0(1), and the call stack may be 0{n) levels deep, so the total space is 0{n). Polynomial 
delay follows as before. As we have seen, the time per child is 0(1), except for the time spent 
sorting the vertices in the fertile_pairs [w] data structure prior to looping over them. That time is 
0{k log k), and we charge the time spent in this step equally to each of the Q{k) children described 
in Lemma|2j In this way, each recursive call gets charged 0(logA:) time, negligible compared to the 
0{k) time the recursive call will spend on listing its own children. □ 

8 Conclusions 

We have provided a general reverse search based framework for generation of all maximal indepen- 
dent sets, and applied it to many important graph classes. Along the way we were led to study new 
dynamic graph data structures for independence and domination. 

One natural problem for additional research is to quantify and reduce the dependence of the 
running time on the sparseness of the graphs in question. Our maximal set generation algorithms 
depend in an exponential way on the sparseness parameter (orientability, degeneracy, or degree) 
of the graphs we consider, due to the way we examine all independent subsets of the sets later (v). 
Can this exponential dependence be reduced? 

Also, to what other non-sparse graph classes can our results be extended? The class of graphs 
that can be ordered so that | later (f)| = 0(1) for all v contains non-sparse graphs, but it is not clear 
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how to perform the dominance or parent tests needed by our algorithms efficiently in such graphs. 
However, we were able to extend our technique to some other graph classes where it did not improve 
previous results. For chordal graphs, if we use an elimination ordering, each independent subset N C 
later (i;) has exactly one vertex, and (by techniques involving counting uniquely dominated vertices) 
we were able to achieve 0(m) time per independent set, matching a prior result of Leung [12]. 
Maximal independent sets for interval graphs may be translated to paths on an associated digraph 
(connect intervals I and J by an edge if J is to the left of I and no interval K is between the 
two), so all maximal independent sets or more general weighted k-hest independent set generation 
problems can be solved in constant time per set by a fc-shortest paths algorithm [7] unrelated to the 
present reverse search approach; a similar approach also works for generating chains and antichains 
in two-dimensional dominance relations. For intersection graphs of disks, balls, squares, or cubes, if 
we order the objects by size then each independent subset of later(i;) must have a constant number 
of objects, so our algorithm can be made to perform in polynomial time per output set, but the 
polynomial appears larger than for general graph independent set enumeration algorithms. 
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